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Abstract:  Upwelling  radiance  measurements  made  with  instruments 
designed  to  float  at  the  sea  surface  are  shaded  both  by  the  instrument 
housing  and  by  the  buoy  that  holds  the  instrument.  The  amount  of 
shading  is  wavelength  dependent  and  is  affected  by  the  local  marine 
and  atmospheric  conditions.  Radiance  measurements  made  with  such 
instruments  should  be  corrected  for  this  self-shading  error  before  be¬ 
ing  applied  to  remote  sensing  calibrations  or  remote  sensing  algorithm 
validation.  Here  we  use  Monte  Carlo  simulations  to  compute  the  self¬ 
shading  error  of  a  commercially  available  buoyed  radiometer  so  that 
measurements  made  with  this  instrument  can  be  improved.  This  ap¬ 
proach  can  be  easily  adapted  to  the  dimensions  of  other  instruments, 
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1  Introduction 

When  an  instrument  is  placed  in  water  to  measure  the  upwelling  light,  the  magnitude 
of  the  local  light  field  is  decreased  by  the  shadow  of  the  instrument.  This  effect  can 
be  large  in  turbid  waters  [1],  causing  several  practical  problems.  Vicarious  calibrations 
of  remote  sensing  radiometers  will  be  in  error  when  performed  with  shaded  in-water 
radiance  measurements,  and  remote-sensing  reflectance  values  will  be  underestimated 
when  calculated  from  shaded  upwelling  radiance  measurements  and  unshaded  down- 
welling  irradiance  measurements.  Furthermore,  because  the  magnitude  of  the  shading 
error  is  wavelength  dependent,  algorithms  that  depend  on  ratios  of  the  radiance  at 
different  wavelengths  to  determine  water  optical  properties  [2],  water  constituent  con¬ 
centrations  [3,  4],  or  bottom  depth  or  features  [5]  will  also  be  in  error.  Therefore,  shad¬ 
ing  corrections  should  be  routinely  applied  to  upwelling  light  measurements.  The  self¬ 
shading  of  in- water  instruments  has  been  investigated  by  Gordon  and  Ding  [1],  Zibordi 
and  Ferrari  [6],  Aas  and  Korsbp  [7],  and  Piskozub,  Weeks,  Schwarz,  and  Robinson  [8]. 

Gordon  and  Ding  [1]  provide  a  semianalytical  model  for  the  self-shading  error  of  sea- 
surface  upwelling  radiance  and  irradiance  measurements  in  optically  deep  waters.  The 
ocean  optics  protocols  for  SeaWiFS  validation  [9]  recommends  the  use  of  this  model,  and 
Zibordi  and  Ferrari  [6]  found  good  agreement  between  the  model  and  their  experimental 
estimates  of  self-shading.  The  fractional  error  in  this  model  is  predicted  with 


£  =  [1  —  exp(— fear)]. 


(1) 
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where  a  is  the  water  beam  absorption  coefficient,  r  is  the  radius  of  the  instrument  hous¬ 
ing,  and  A:  is  a  constant  that  depends  on  the  measurement  type  and  on  the  illumination 
conditions.  Gordon  and  Ding  [1]  derived  values  of  k  with  Monte  Carlo  simulations  of  a 
two-dimensional  shading  disk  placed  just  below  the  sea  surface.  For  radiance  measure¬ 
ments  and  a  sun  in  a  black  sky,  k  ~  2/tan0ou)i  where  9qw  is  the  in- water  solar  zenith 
angle.  The  shading  error  for  general  sky  conditions  can  be  approximated  with  [1,  6] 


£  = 


Sskyf 


1  +  / 


(2) 


where  £sun  is  the  shading  error  for  direct  sunlight  from  the  appropriate  sun  angle,  £sky  is 
the  shading  error  for  skylight,  and  /  is  ratio  of  the  downwelling  irradiance  from  skylight 
to  that  from  direct  sunlight. 

Unfortunately,  for  upwelling  light  sensors  mounted  in  instruments  that  use  a  buoy 
to  float  on  the  sea  surface,  Eq.  (1)  is  inappropriate.  Because  the  buoy  has  a  larger 
radius  than  the  sensor  housing  and  is  vertically  separated  from  the  sensor  depth,  it 
could  be  expected  that  the  actual  shading  of  the  instrument  in  optically  deep  water  is 
greater  than  that  predicted  by  Eq.  (1)  for  the  sensor  housing  alone  but  less  than  that 
predicted  by  Eq.  (1)  for  the  buoy.  Equation  (1)  is  also  not  valid  for  optically  shallow 
waters,  in  which  the  shading  error  is  highly  dependent  on  water  depth  and  the  seafloor 
characteristics. 

Here  we  use  Monte  Carlo  simulations  to  quantify  the  shading  effects  of  cylindrical 
instruments  and  the  additional  shading  of  a  wider  cylindrical  surface  buoy.  Our  specific 
objective  is  to  derive  a  shading  correction  algorithm  for  the  Hyper-TSRB  (Satlantic 
Inc.,  Halifax,  Nova  Scotia,  Canada),  which  was  extensively  used  in  recent  CoBOP  [10] 
experiments  for  in  situ  studies  and  as  sea-truth  for  aircraft  remote  sensing  systems.  The 
approach  demonstrated  here  and  most  of  our  qualitative  conclusions  apply  to  any  set 
of  instrument  dimensions;  only  the  quantitative  results  are  specific  to  the  TSRB. 


2  Methods 

Simulated  responses  of  an  upwelling  radiance  detector  with  and  without  the  presence 
of  a  shading  object  were  computed  using  Backward  Monte  Carlo  (BMC)  ray  tracing 
as  described  in  Ref.  11.  In  BMC  codes  simulated  photons  are  emitted  from  the  sensor 
and  traced  backwards  to  the  light  source.  For  geometries  involving  an  extended  source 
(the  sky  in  our  case)  and  small  sensor  (as  on  the  TSRB),  the  BMC  technique  is  numer¬ 
ically  efficient  because  every  traced  photon  contributes  to  the  statistical  estimate  of  the 
upwelling  radiance  at  the  instrument.  The  radiance  sensor  is  treated  as  a  point  sensor 
located  at  the  bottom  center  of  a  cylindrical  casing.  After  emission  from  the  sensor, 
the  photons  are  traced  using  standard  techniques  for  forward  Monte  Carlo  ocean  optics 
simulations  as  described  in  Ref.  12.  Shading  is  simulated  by  removing  any  photon  whose 
path  intersects  the  cylinder.  The  equations  necessary  to  check  for  photon-cylinder  in¬ 
tersections  in  a  Monte  Carlo  algorithm  are  given  later  in  this  section.  In  our  BMC  code 
the  effects  of  the  illumination  conditions  are  applied  last,  allowing  us  to  apply  many 
different  sky  conditions  to  the  same  in- water  ray-tracing  results.  For  most  computations 
we  used  the  Fournier-Forand  [13]  scattering  phase  function  with  the  parameters  chosen 
to  closely  match  the  Petzold  particle  [14]  scattering  phase  function.  The  seafloor  was 
assumed  to  be  Lambertian  with  a  given  irradiance  reflectance,  and  the  shape  of  the  sea 
surface  was  taken  to  be  that  of  capillary  waves  as  described  by  Cox-Munk  statistics  [15]. 

The  TSRB  was  modeled  as  two  cylinders  extending  into  the  water  from  the  sea 
surface.  The  main  body  has  a  4.4  cm  radius  and  extends  66  cm  into  the  water,  whereas 
the  buoy  has  a  15  cm  radius  and  extends  12  cm  into  the  water.  Potential  shading  due 
to  the  portion  of  the  TSRB  that  is  above  the  waterline  was  not  included,  but  we  expect 
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this  to  be  insignificant  for  upwelling  radiance  measurements.  The  TSRB  radiance  sensor 
response  characteristics  were  modeled  from  measurements  provided  by  Satlantic  [16]. 

The  shading  error  s  was  computed  for  each  set  of  water  properties  and  illumination 
conditions  with 

£  =  (3) 

where  is  the  simulated  sensor  response  with  shading  and  is  the  simulated 

unshaded  sensor  response.  Many  runs  of  the  Monte  Carlo  code  were  made  for  each 
case  but  with  different  seeds  of  the  random  number  generator  so  that  the  statistical 
uncertainty  [17]  of  the  simulations  could  be  estimated. 

The  parametric  equation  for  the  line  segment  describing  a  photon  path  is 


X 

Xl 

a 

y 

= 

yi 

-p 

p 

z 

■^1 

.  O'  . 

(4) 


where  x  and  y  are  horizontal  positions,  z  is  depth  (m),  (xi,  yi,  zi)  is  the  starting  point 
of  the  line  segment,  a,  j3,  and  7  are  the  direction  cosines  of  the  line  segment,  and  L  is 
the  length  of  the  line  segment  (i.e.,  the  distance  the  photon  travels).  The  parametric 
equation  for  the  solid  cylinder  describing  the  instrument  is 


X 

Xo 

pcos{9) 

Q  <  p  <r 

y 

= 

yo 

+ 

psm{0) 

,  0  <  6»  <  27r, 

(5) 

z 

Zo 

u 

0  <  u  <  h, 

where  (xq,  yo,  zq)  is  the  location  of  the  bottom  center  of  the  cylinder  and  r  and  h  are 
the  cylinder  radius  and  height.  The  intersection  points  between  the  line  coincident  with 
the  photon  path  and  an  infinite  cylinder  coincident  with  the  cylinder  wall  satisfy 


(xi  -  xo  +  asf  +  (j/i  -  2/0  +  l3sY  =  r'^- 


(6) 

(7) 


Equation  (6)  is  quadratic  in  s,  the  solutions  of  which  are 


s  = 


-B  ±  -  AAC 

2^4 


A  = 

B  =  2a(xl  —  xO)  +  2/3(yl  —  yO), 

C  =  (xl  —  xO)^  +  (yl  —  yO)^  —  r^. 


(8) 


There  are  no  intersection  points  if  {B^  —  4ylC)  <  0.  Otherwise,  for  solutions  si  and 
S2,  the  range  of  values  of  s  for  which  the  infinite  line  intersects  an  infinitely  tall  solid 
cylinder  is  si  <  s  <  S2.  The  finite  photon  path  (of  length  L)  is  a  subsegment  of  this 
infinite  line.  Values  of  s  that  satisfy  s  €  [max(0,  s)  <  s  <  min(L,  S2)]  define  the  portion 
of  the  photon  path  that  intersects  the  infinitely  tall  solid  cylinder.  If  this  range  of  s  is 
finite,  we  find  the  boundaries  of  the  corresponding  range  of  u  by  substituting  the  limits 
of  s  into  [from  Eqs.  (4)  and  (5)] 


u  =  zi  -  zo  +  7s.  (9) 

If  this  range  of  u  overlaps  0  <  u  <  h  at  any  point,  then  the  photon  is  known  to  hit  the 
cylinder  of  height  h. 

3  Results 

Figure  1  shows  the  TSRB  self-shading  error  in  optically  deep  water  versus  the  water 
beam  absorption  coefficient  a  (m“^)  compared  with  the  errors  predicted  by  Eqs.  (1) 
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Fig.  1.  Self-shading  error  of  the  TSRB  (solid  line)  in  deep  water  for  hja  =  2  com¬ 
pared  with  that  computed  with  Eqs.  (1)  and  (2)  for  radii  0.044  m  and  0.15  m.  The 
TSRB  error  was  computed  for  an  empirical  sky  model  at  480  nm  with  the  sun  zenith 
angle  9q  =  30  degrees. 

and  (2)  for  shading  disks  with  radii  of  the  TSRB  instrument  body  and  TSRB  buoy, 
respectively.  As  expected,  the  shading  error  of  the  TSRB  is  larger  than  that  of  a  disk 
with  r  =  0.044  m  (the  TSRB  body  radius)  because  of  the  presence  of  the  buoy,  but 
is  less  than  that  of  a  disk  with  r  =  0.15  m  (the  buoy  radius)  because  there  is  vertical 
separation  between  the  buoy  and  the  sensor.  Because  the  optical  distance  between  the 
buoy  and  the  sensor  is  larger  in  turbid  waters  than  in  clear  waters,  the  TSRB  shading 
error  in  turbid  water  is  closer  to  that  of  a  disk  with  the  instrument  radius  than  that  of 
a  disk  with  the  buoy  radius  and  vice  versa  in  clear  waters. 

Figure  2  shows  the  shading  error  versus  a  for  optically  deep  water,  b/a  =  2,  and 
five  solar  zenith  angles  Oq  for  a  sun  in  a  black  sky.  For  large  values  of  9o  the  error 
increases  approximately  linearly  with  a  for  0.02  <  a  <  1,  whereas  for  smaller  sun  angle, 
and  therefore  larger  shading  errors,  the  error  depends  on  a  in  a  nonlinear  manner. 
Equation  (1)  gives  a  poor  fit  to  the  curves  in  Fig.  2,  either  underestimating  the  error  at 
low  values  of  a  or  overestimating  the  error  at  high  values  of  a,  depending  on  the  chosen 
value  of  k. 

As  shown  in  Fig.  2,  self-shading  error  depends  strongly  on  sun  position.  The  top  plot 
in  Fig.  3  shows  the  upwelling  radiance  L„  versus  solar  zenith  angle  both  when  the  TSRB 
is  present  and  in  its  absence.  These  simulations  are  for  a  sun  in  a  black  sky,  optically 
deep  water,  and  a  =  b  —  0.2  m“^.  The  values  of  shown  in  Fig.  3  are  normalized 
by  the  downward  irradiance  at  the  sea  surface  Eci(0).  The  shading  error,  shown  in  the 
bottom  plot  of  Fig.  3,  is  large  for  small  sun  angles  and  small  for  large  sun  angles.  Here, 
too,  we  see  that  the  TSRB  shading  error  in  deep  water  lies  between  those  predicted 
by  Eq.  (1)  for  the  two  TSRB  radii.  Also  shown  in  Fig.  3  is  the  error  due  to  the  main 
body  of  the  TSRB  without  the  buoy.  It  can  be  seen  that  the  buoy  has  a  large  effect,  but 
only  for  9q  <  20°.  This  was  found  to  be  true  in  general,  and  therefore  the  shading  error 
of  the  TSRB  has  a  much  stronger  dependence  on  Oq  than  does  that  of  a  non-buoyed 
cylindrical  instrument.  Figure  4  shows  the  TSRB  shading  error  versus  solar  zenith  angle 
for  a  sun  in  a  black  sky,  six  values  of  a,  and  b  —  2a.  For  all  values  of  a,  the  shading  error 
is  significant  at  small  values  of  9o  and  decreases  with  increasing  Oq. 

For  a  given  value  of  a  the  shading  error  in  optically  deep  water  is  generally  smallest 
for  large  values  of  b  because  scattering  enables  the  regions  shadowed  from  direct  light 
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Fig.  2.  Self-shading  error  of  the  TSRB  in  deep  water  for  hja  =  2,  five  solar  zenith 
angles,  and  no  skylight. 


X  10'^ 


Fig.  3.  Shaded  and  unshaded  normalized  radiance  and  percent  shading  error  of  the 
TSRB  and  of  the  TSRB  cylindrical  body  without  the  buoy  compared  with  Eq.  (1) 
for  radii  0.044  m  and  0.15  m.  The  computations  are  for  a  sun  in  a  black  sky,  deep 
water,  and  a  =  b  =  0.2  m~^. 
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Fig.  4.  Percent  shading  error  of  the  TSRB  in  deep  water  versus  solar  zenith  angle 
Bq  for  a  sun  in  a  black  sky,  six  values  of  a,  and  bja  =  2. 


Fig.  5.  Percent  shading  error  of  the  TSRB  in  deep  water  for  a  =0.02  (dashed), 

0.05  (solid),  and  1.0  (dotted)  and  9o  =  0°  (top,  blue),  10°  (middle,  green), 
and  20°  (bottom,  red). 

to  be  filled  in  with  scattered  light.  However,  we  found  this  effect  to  be  significant  only 
for  small  values  of  Oq,  where  the  shading  effect  is  greatest.  Shown  in  Fig.  5  is  the  TSRB 
shading  error  versus  b  for  Oq  =  0,  10,  and  20  degrees,  each  for  a  =  0.2,  0.5,  and  1.0  m“^. 
It  can  be  seen  that  the  slope  of  shading  error  curves  are  similar  for  constant  9q  but  that 
this  slope  decreases  with  increasing  Oq. 

Table  1  lists  the  TSRB  shading  error  for  flat,  optically  deep  water  both  for  a  diffuse 
sky  and  for  a  sun  in  a  black  sky  at  eight  solar  zenith  angles.  Results  are  shown  for  sbc 
values  of  a  and,  for  6q  <  20°,  three  values  of  b/a.  For  6q  >  20°  the  dependence  on  b  was 
found  to  be  insignificant  and  the  results  shown  in  Table  1  are  those  for  bja  =  1,  2,  and 
4  averaged  together.  Estimates  of  the  statistical  uncertainty  in  the  error  values,  derived 
from  the  standard  deviation  of  the  simulations,  is  provided  for  the  smallest  and  largest 
values  of  a. 

Wind-induced  surface  capillary  waves  tend  to  spread  out  the  angular  distribution  of 
incoming  light.  This  can  reduce  the  shading  error  when  the  sun  is  high  in  the  sky  by 
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Table  1.  Percent  shading  error  (lOOxe)  of  a  TSRB  for  given  values  of  absorption 
coefficient  a,  scattering  coefficient  b,  and  solar  zenith  angle  Oq. 


do 

a  = 

0 

02  m"^ 

0.05 

0.1 

0.2 

0.5 

1 

.0 

0° 

bja- 

=  1 

6.1 

± 

0.6% 

11.7 

22.0 

35.8 

55.5 

69.5 

± 

1.4 

2 

6.0 

± 

0.7 

11.1 

20.8 

33.3 

49.3 

58.1 

± 

1.0 

4 

5.9 

± 

0.8 

11.0 

18.0 

26.3 

37.3 

43.0 

± 

0.5 

10° 

bja  = 

=  1 

3.2 

± 

0.3 

7.6 

13.1 

21.6 

39.9 

53.6 

± 

0.5 

2 

3.2 

± 

0.2 

6.9 

12.3 

20.4 

35.0 

46.5 

± 

0.4 

4 

2.6 

± 

0.3 

6.4 

11.0 

17.7 

28.6 

36.8 

± 

0.3 

o 

O 

b/a  = 

=  1 

0.6 

± 

0.2 

2.0 

4.4 

8.2 

17.4 

29.8 

± 

0.4 

2 

0.9 

± 

0.2 

2.2 

4.1 

7.6 

17.3 

28.8 

± 

0.4 

4 

1.0 

± 

0.3 

2.0 

4.3 

7.7 

16.4 

26.3 

± 

0.2 

o 

O 

CO 

b/a- 

=  2 

0.5 

± 

0.1 

1.2 

2.2 

4.9 

11.3 

20.9 

± 

0.2 

o 

o 

bja  = 

=  2 

0.5 

± 

0.1 

0.7 

1.8 

3.5 

8.8 

16.4 

± 

0.1 

50° 

b/a  - 

=  2 

0.1 

± 

0.1 

0.6 

1.1 

2.8 

6.9 

13.9 

± 

0.1 

o 

O 

b/a  - 

=  2 

0.2 

± 

0.1 

0.6 

1.0 

2.4 

6.1 

12.1 

± 

0.1 

70° 

b/a  - 

=  2 

0.0 

± 

0.1 

0.7 

1.1 

2.3 

5.3 

11.2 

± 

0.2 

diffuse 

b/a- 

=  1 

0.5 

± 

0.1 

1.3 

2.6 

5.1 

11.4 

20.0 

± 

0.2 

2 

0.6 

± 

0.1 

1.4 

2.5 

4.9 

10.9 

19.4 

± 

0.2 

4 

0.5 

± 

0.1 

1.5 

2.4 

4.7 

10.2 

18.4 

± 

0.1 

bending  nearby  incident  light  in  toward  the  sensor.  However,  we  found  this  effect  to  be 
insignificantly  small.  Apparently,  some  photons  are  refracted  toward  the  shaded  region 
and  some  are  refracted  away,  giving  little  overall  effect. 

In  general  the  upwelling  radiance  measurement  is  composed  of  light  that  is  scattered 
by  the  water  column  into  the  sensor  field  of  view  and  of  light  that  is  reflected  off  the 
seafloor.  The  former  dominates  in  deep  waters  with  highly  absorbing  seafloors  and  the 
latter  dominates  in  shallow  waters  with  highly  reflective  bottoms.  The  shading  error  in 
the  water  column  component  is  approximately  that  given  in  Table  1  for  large  depths,  but 
increases  as  the  water  depth  is  decreased.  The  component  of  the  radiance  measurement 
due  to  reflectance  off  the  bottom  also  contains  some  shading  error  because  the  TSRB 
casts  a  shadow  on  the  bottom  that  will  always  be  at  least  partially  in  the  field  of  view 
of  the  radiance  sensor.  However,  we  found  that  usually,  but  not  always,  the  shading 
error  of  the  bottom  component  is  smaller  than  that  for  the  water  column.  Therefore, 
for  shallow  waters  and  highly  reflecting  bottoms,  where  the  bottom  signal  dominates 
over  the  water  column  signal,  the  overall  shading  error  is  often  much  smaller  than  in 
optically  deep  water  with  the  same  water  properties.  For  example.  Fig.  6  shows  the 
shading  error  as  a  function  of  total  water  depth  for  the  case  of  a  =  0.2  m“^,  b  =  0.4 
m“^,  and  bottom  irradiance  reflectance  Rb  =  0.2.  In  very  shallow  water  the  radiance 
measurement  is  dominated  by  the  component  of  light  being  reflected  off  the  bottom. 
As  the  depth  is  increased,  the  shading  error  in  the  radiance  reflected  off  the  bottom 
decreases,  as  does  the  overall  shading  error.  However,  as  the  water  depth  is  further 
increased  beyond  4  m,  the  water-column  component  of  the  upwelling  radiance  begins  to 
dominate  and  the  overall  shading  error  increases  toward  the  optically  deep  values. 

4  Discussion 

The  self-shading  of  in-water  light-measuring  instruments  depends  on  the  illumination 
conditions,  water  depth,  and  water  and  seafloor  optical  properties.  The  shading  of  sea- 
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Fig.  6.  TSRB  self-shading  error  as  a  function  of  water  depth  for  absorption  coeffi¬ 
cient  a  =  0.2  m“^,  scattering  coefficient  b  =  0.4  and  bottom  albedo  Rf,  =  0.2. 

The  dashed  lines  show  the  shading  error  in  optically  deep  waters  for  solar  zenith 
angles  =  0°,  10°,  and  20°. 


surface  upwelling  radiance  and  irradiance  measurements  made  with  cylindrical  instru¬ 
ments  in  optically  deep  water  can  be  estimated  with  Eq.  (1).  However,  for  other  config¬ 
urations  and  environmental  conditions  an  empirical  approach  is  required.  From  three- 
dimensional  Monte  Carlo  simulations  we  found  that  the  total  shading  error  of  a  buoyed 
instrument  in  optically-deep  water  lies  between  that  predicted  by  Eq.  (1)  for  the  radius 
of  the  instrument  body  and  that  predicted  by  Eq.  (1)  for  the  buoy  radius.  The  error  is 
closer  to  the  former  in  turbid  waters  and  closer  to  the  latter  in  clear  waters.  Because  the 
shading  error  is  greatest  at  small  solar  zenith  angles,  upwelling  radiance  measurements 
should  be  taken  when  the  solar  zenith  angle  is  greater  than  20°,  if  possible. 

Upwelling  radiance  measurements  obtained  with  a  buoyed  instrument  should  be 
corrected  with  data  similar  to  that  presented  in  Table  1.  The  following  steps  can  be 
used  to  apply  this  correction: 


1. 


Determine  the  value  of  a(A)  for  each  radiance  wavelength  A.  Ideally  a(A)  is  meas¬ 
ured  from  water  samples  [19,  20]  or  with  an  ac-9  (WET  Labs,  Inc.,  Philomath, 
Oregon).  Otherwise,  given  upward  and  downward  irradiance  measurements  i?„(A) 
and  Ed{X)  at  two  or  more  depths  near  the  surface,  a(A)  can  be  estimated  with, 
for  example  [21], 


jLrf(A)[l-i?(A)]cosgo 
0.6-b  [0.47-b2.5i?(A)]cos6»o’ 


(10) 


where  R  =  Eu/Ed  and  Kd  =-d(log(ifd))/dz.  Alternately,  given  values  of  L„(A)  at 
more  than  one  depth,  a  can  be  estimated  with  [7]  a(A)  =  x^lW,  where  =- 
d(log(L„))/dz  and  y  is  a  constant  of  proportionality  valid  for  the  local  region. 
Given  a(A)  at  several  wavelengths,  a(A)  can  be  modeled  for  all  wavelengths  with 
a  bio-optical  model  [12].  In  the  absence  of  any  ancillary  data,  use  an  a(A)  that  is 
at  least  as  large  as  the  published  value  for  pure  seawater  [12,  22]. 


2.  Estimate  6(A)  for  each  wavelength.  This  is  only  important  if  a(A)  is  large  and  6o 
is  small.  Otherwise,  let  6(A)  =  2a(A).  Pure  water  values  of  6  are  given  in  Refs.  12, 
22,  and  23,  and  bio-optical  models  for  6(A)  are  given  in  Ref.  12. 


3.  Given  a(A)  and  6(A),  interpolate  Table  1  to  obtain  esun(A)  for  the  appropriate  sun 
zenith  angle  Oq.  Also  interpolate  the  “diffuse”  portion  of  Table  1  to  obtain  esky(A). 


#32933  -  $15.00  US 
(C)  2001  OS  A 


Received  April  17,  2001;  Revised  May  01,  2001 
7  May  2001  /  Vol.  8,  No.  10  /  OPTICS  EXPRESS  569 


Fig.  7.  Example  TSRB  shading  correction.  The  upper  plot  shows  the  water  absorp¬ 
tion  spectrum  a(A),  the  ratio  /(A)  of  skylight  to  sunlight,  and  the  corresponding 
shading  error  £(A).  The  lower  plot  shows  the  measured  and  corrected  upwelling 
radiance  spectra. 

For  each  wavelength  determine  the  ratio  of  skylight  to  direct  sunlight,  /(A).  This 
can  be  done  from  direct  measurement  (e.g.,  with  a  sun  photometer),  with  an 
empirical  atmospheric  model  [18,  24],  with  the  procedure  given  by  Gordon  [25], 
or  by  selecting  published  values  [26]. 

Determine  the  total  shading  error  e(A)  for  each  wavelength  with  Eq.  (2). 

Compute  the  true  radiance  spectrum  from  the  measured  radiance  T™(A) 

with 

L^nX)  =  L^{X)/{l-e{X)).  (11) 

An  example  TSRB  shading  correction  is  demonstrated  in  Fig.  7.  The  water  absorp¬ 
tion  spectrum  shown  in  the  upper  plot  is  typical  of  case  I  waters  with  relatively  high 
chlorophyll  concentration  (10  mg/m)  [12].  The  measured  radiance  spectrum  shown  in 
the  lower  plot  is  for  a  solar  zenith  angle  of  10°  and  a  sky  with  20%  clouds.  The  ratio 
of  skylight  to  sunlight  /(A)  shown  in  the  upper  plot  was  computed  using  an  empiri¬ 
cal  atmospheric  model  [24].  The  corresponding  shading  error  e(A)  was  determined  with 
Table  1  and  Eq.  (2)  and  was  applied  to  the  measured  radiance  spectrum  to  obtain  the 
corrected  radiance  spectrum  shown  in  the  lower  plot. 

A  Fortran  routine  that  replaces  steps  3-5  above  is  available  from  the  authors.  This 
routine  uses  neither  Table  1  nor  Eq.  (2),  but  instead  applies  modeled  sky  conditions 
to  tabulated  values  of  BMC  ray-tracing  statistics.  It  is  important  to  note  that  the 
shading  error  values  given  in  Table  1  are  only  valid  for  optically  deep  water,  calm  sea- 
surface  conditions,  and  particular  TSRB  dimensions.  Different  tables  need  to  be  used 
for  different  situations.  Although  it  is  not  practical  to  provide  many  tables  here,  the 
shading  error  for  many  different  water  depths,  bottom  types,  and  sea  surface  roughness 
are  being  compiled  for  use  in  the  Fortran  routine. 

Acknowledgments 

This  work  was  supported  by  the  U.S.  Office  of  Naval  Research.  Author  C.D.M.  was 
supported  by  ONR’s  Environmental  Optics  Program.  We  thank  Marcos  Montes  for  his 
helpful  suggestions. 


4. 

5. 

6. 


#32933  -  $15.00  US 
(C)  2001  OS  A 


Received  April  17,  2001;  Revised  May  01,  2001 
7  May  2001  /  Vol.  8,  No.  10  /  OPTICS  EXPRESS  570 


